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1972 ANNUAL REPORT ON THE INVESTIGATION OF 
CRITICAL PRESSURE BURNING OF FUEL DROPLETS 


Summary 


This report discusses activities under NASA Contract NGR 39-009-077 
for the period January 1, 1972, through December 31, 1972. During this 
period the investigation was divided into four phases, the results for 
each phase may be summarized as follows: 

1. Droplet Combustion at High Pressures . Combustion of methanol, 
ethanol, propanol- 1, n-pentane, n-heptane and n-decane from porous 
spheres was observed in air under natural convection conditions at 
pressures up to 100 atm. The pressure levels of the tests were 
high enough so that near critical combustion was observed for 
methanol and ethanol. Measurements were made of the burning rate 
and liquid surface temperatures of the fuels. Due to the high 
pressures, the phase equilibrium models of the analysis included 
both low and high pressure versions, the latter allowing for real 
gas effects and the solubility of the combustion product gases 

in the liquid phase. The burning rate predictions of the various 
theories were similar and in fair agreement with the data. The 
experiments indicated the approach of critical burning conditions 
for methanol and ethanol at pressures on the order of 80-100 atm, 
which was in good agreement with the predictions of both the low 
and high pressure analysis for these fuels. Critical burning 
conditions could not be approached for the remaining fuels due to 
the formation of soot deposits on the sphere at pressures in the 
range 30-60 atm. 

2. High Pressure Droplet Burning in Combustion Gases . An experimental 
apparatus was constructed to allow measurements of the steady com- 
bustion of liquid fuels from porous spheres in an environment that 
simulates combustion chamber conditions. The apparatus employs a 
burner contained within a pressure vessel to allow operation at 
various oxygen concentrations, gas velocities and pressures 
ranging up to 35 atm. During the coming report period, testing 
and analysis will be completed for the same range of fuel types 
considered in the natural convection tests. 

3. Liquid Strand Combustion . Burning rates and liquid surface temp- 
eratures were measured for burning liquid strands of hydrazine over 
the pressure range 7-600 psia. At pressures greater than 1 atm, 
the burning rate was found to increase linearly with pressure. The 
addition of water to the hydrazine resulted in a marked reduction 
in the burning rate, Liquid surface temperatures were found to be 
in good agreement with predictions based on a low pressure phase 
equilibrium model for pressures up to 300 psia. 
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4. Oscillatory Combustion . An experimental apparatus has been con- 
structed to allow measurement of the response of a burning liquid 
hydrazine strand to imposed pressure oscillations. Due to the ability 
of hydrazine to burn at very slow rates as a strand, this apparatus 
is capable of observing the response at frequencies comparable to 
the critical frequency of liquid phase transient effects, providing 
a good test of the oscillatory theory. Equations are developed for 
the oscillatory theory of the response through first order in the 
amplitude of the pressure oscillations. During the coming report 
period experimental results will be completed and compared with 
theoretical computations. 



I . Introduction 


The objective of this study is to investigate the combustion and 
evaporation of liquid fuels at high pressures. Particular emphasis is 
placed on conditions where the liquid surface approaches the thermodynamic 
critical point during combustion. The influence of transient effects on a 
burning liquid fuel is also being investigated through both analysis and 
measurements of the response of liquid monopropellant combustion to imposed 
pressure oscillations. This report gives a summary of progress on the 
investigation for the period January 1, 1972, to December 31, 1972. 

Work during this report period was divided into four phases as follows: 

(1) Droplet Combustion at High Pressures, which considers both measurement 
and analysis of the porous sphere burning rate of liquids in a natural con- 
vection environment at elevated pressure. Work in this area was concluded 
during the report period. (2) High Pressure Droplet Burning in Combustion 
Gases, which involves steady burning and evaporation of liquids from porous 
spheres in a high pressure environment that simulates actual combustion 
chamber conditions. Work in this area is currently in progress. (3) Liquid 
Strand Combustion, which considers the burning rate, the state of the liquid 
surface and the liquid phase temperature distribution of a burning liquid 
monopropellant column over a range of pressures. Work in this area was con- 
cluded during the report period. (4) Oscillatory Combustion, which is a 
theoretical and experimental investigation of the response of a burning 
liquid monopropellant to pressure oscillations. Work in this area is currently 
in progress. 

In the following sections, the activity in each phase is summarized 
for the report period. Wherever possible, the details of the work are 
presented by reference to past publications and only material not pre- 
viously reported is given extended coverage. Reference 1, 15, and 25 
comprise other publications resulting from work done on this project which 
have appeared during the period of this report. 


II . Droplet Combustion at High Pressures 


Introduction: 

The objective of this part of the investigation was to study high 
pressure droplet combustion for a variety of fuels. This included the 
experimental determination of burning rates and liquid surface tempera- 
tures. The experimental results were compared with droplet combustion 
theories which both neglected and considered real gas effects. The fuels 
used in the tests consisted of methanol, ethanol, propanol-1, n-pentane, 
n-heptane, and n-decane= The pressure range varied from 1 to 51 atmospheres 
for the hydrocarbon fuels and from 1 to 78 atmospheres for the alcohols. 
Porous alundum spheres were used to simulate the fuel droplets. The sphere 
sizes included 0.63 cm. 0.95 cm and 1.9 cm diameter spheres. 

In addition to the determination of burning rates and liquid surface 
temperatures, a series of dark field photographs were taken. Methanol and 
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ethanol were the fuels used in these tests. The sphere sizes employed were 
0,63 cm, 0.95 cm and 1,9 cm, in the pressure range from 1 to 35 atmospheres. 

The findings of this portion of the investigation are summarized in greater 
detail in Reference 1. 

Experimental Apparatus: 

A sketch of the overall system is shown in Figure 1. The high-pressure 
combustion apparatus consists of a thick-walled steel cylindrical vessel 
66 cm, long with an inside diameter of 13 cm. The chamber is fitted with 
windows to allow observation of the combustion process. 

The porous spheres used in the combustion tests are made of alundum. 

The fuel is fed to the center of the sphere through a stainless steel 
water cooled tube. The outside diameter of the coolant tube was 0,20 cm 
for the 0,64 cm diameter sphere and 0.32 cm for the larger spheres. The 
fuel flow rate was controlled by a variable displacement diaphragm pump 
and measured with a system of graduated burets at the pump inlet. The 
steady burning rate was determined as the flow rate where the surface of 
the sphere is fully wet but not dripping. 

Air for the combustion process is supplied by a 3000 psia compressor 
which is accurately controlled by a critical flow jeweled orifice system 
of the type described by Anderson and Friedman ( 2 ) . The hot exhaust gases 
leaving the reactor are cooled with water in a concentric tube heat ex- 
changer. Water condensate is collected in a trap which is periodically 
blown off to a drain. 

Liquid surface temperatures were measured using two 0.0076 cm diameter 
chromel-alumel thermocouples. All liquid surface temperature determinations 
were made using the 0.95 cm diameter sphere. The thermocouples were cemented 
flush with the surface and located 60 degrees apart along the periphery of 
the sphere, A constant gas composition was maintained around the sphere by 
admitting air through a manifold at the bottom of test chamber. The drift 
velocity of the flow of air past the sphere was sufficiently low so that 
natural convection was the predominant convection effect. 

The dark field photographs were taken using a 4 in. x 5 in. Super Graflex 
camera fitted with a 135 mm Optar lens. The film used was Polaroid type 57 
land film with an ASA rating of 3000, 

Theoretical Analysis: 

The present theory is similar in many respects to that of References 3 and 
4 for high pressure combustion. The major point of difference involves the 
different boundary conditions at the liquid surface for the present porous 
sphere combustion case as opposed to steady droplet combustion, and a 
correction to the burning rate for natural convection effects. 

The theory is divided into a gas phase model of the combustion process 
and a phase equilibrium model for conditions at the liquid surface. In the 
gas phase model, the effect of convection is treated by the use of the 
multiplicative correction of the burning rate predicted in the absence of 
convection similar to that employed in References 5 and 6, Therefore, the 
basic analysis assumes spherical symmetry and neglects convection effects. 

The process is assumed to be steady and the total pressure is taken to be 
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uniform throughout the system. Both dissociation and radiation are neglected 
and the reaction is assumed to be confined to an infinitely thin flame surface 
where fuel and oxidizer combine in stoichiometric proportions. 

Only concentration diffusion is considered in the gas phase and the 
influence of compressibility on transport properties is neglected. In 
addition, the concentration dependence of the thermal conductivity is 
neglected and the binary diffusivities of all species are taken to be 
the same, although different values of each of these properties can be 
employed inside and outside the flame. The common unity Lewis number 
assumption was not employed in the analysis since earlier studies (3, 4) 
have shown that the value of the Lewis number has a strong influence on 
conditions at the liquid surface. 

A sketch of the model of the burning droplet is shown in Figure 2. 

In the present experiment, the liquid fuel was pumped from a storage vessel 
at atmospheric pressure to the center of the sphere. Therefore, due to 
the low solubility of gases in the test fuels at low pressures, it is 
appropriate to assume that the liquid entering the sphere has a neglibible 
dissolved gas concentration. Under this assumption, the liquid phase flux 
of dissolved gas is zero and the fuel is the only component with a finite 
molar flux inside the flame surface. 

The specific heat of each species was assumed to be a linear function 
of temperature* 


C . 
Pi 


a. + b. T „ 

l l 


CD 


where a^ and b^ are fixed constants for species i. 

The thermal conductivity was also assumed to be proportional to 
temperature in the regions inside and outside the flame 

.A = A s (T/T s ). 


( 2 ) 


The quantity 


A/(cD) 


( 3 ) 


is only a weak function of temperature and composition and was assumed 
to be a constant inside and outside the flame. 


With due allowance for the fact that the fuel mole flux fraction 
is unity in the region inside the flame, the details of the gas phase 
analysis are very similar to that presented in References 3 and 4. The 
specific results obtained in the present case are given in Reference 1. 
In addition to this variable property - variable specific heat analysis, 
calculations were also performed for variable property - constant 
specific heat and constant property models. 


*See the nomenclature at the end of this section for the meaning of symbols. 
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The semi-empirical corrections for the effect of natural convection 
on the burning rate, 6 was based on the natural convection heat transfer 
correlation determined by Yuge. 7 In addition, the Grashoff number for 
burning spheres suggested by Spalding ® was employed in the correlation. 

The specific equation used in the calculations is as follows 


ft i n nm n 1/3, ,3 / 2,1/4 

-r— = 1 + 0.221 Pr (d g/v ) (4) 

o 

In order to avoid ambiguity, the properties used in this correlation were 
taken at the ambient conditions of the burning sphere. 

Three models were employed for computing phase equilibrium at the 
liquid surface. The simplest model (low pressure theory) neglected the 
solubility of the product gases in the liquid phase and other high pres- 
sure corrections. The fuel mole fraction at the' liquid surface was taken 
to be the vapor pressure of the pure fuel, at the surface temperature, divided 
by the total pressure. The total enthalpy rise of vaporization was determined 
by summing the compressed liquid enthalpy change at T (the fuel inlet 
temperature), the heat of vaporization at T , the ideal gas enthalpy rise 
between T and T (the surface temperature) and the enthalpy deviations 
between saturated vapor and the ideal gas states at T and T _ . In this 
case, the ideal gas enthalpy rise was computed by integrating actual specific 
heat correlations between T and T , as opposed to the linearized specific 
heat correlation (Equation °1) employed in the gas phase analysis. The 
enthalpy deviations were obtained from the tables of Lyderson, et al, pre- 
sented in Reference 9. 

The high-pressure theories considered solubility and other high-pres- 
sure effects through the use of a modified Redlick-Kwong equation of state 
with mixing rules given by Prausnitz and Chueh. u For combustion in air, 
the major gaseous species at the liquid surface are fuel vapor, nitrogen, 
carbon dioxide and water vapor. Since nitrogen predominates the non-fuel 
gases in this system, a simplified version of the high pressure theory as- 
sumed that this system could be represented by a binary mixture of fuel and 
nitrogen. The more complete version of the high pressure theory considered 
the complete quaternary system; fuel, nitrogen, carbon dioxide and water. 

The details of the formulation of the high pressure phase equilibrium 
model are given in Reference 4. The dimensionless constants required by 
the Redlick-Kwong equation of state were obtained by setting the' first and 
second isothermal derivatives of pressure, with respect to volume, equal to 
zero at the critical point. The binary interaction parameters, k.,, required 
by the theory are listed in Reference 3 for the paraffins and the 1 iombustion 
product gases. For the alcohols, the k. . values were taken to be the same 
as that of the hydrocarbon homomorph of 1 iuel. These values are listed in 
Reference 1, For the high pressure theories, the enthalpy deviations required 
in the calculation of the total heat of vaporization of the fuel were computed 
directly from the Redlick-Kwong equation of state. The remaining terms com- 
prising the enthalpy rise of vaporization were computed in the same manner as 
in the low pressure theory calculations. 

The sources of property data and correlations were the same as in earlier 
studies, 3, 4, 5 xh e particular values employed in the present calculations 
are given in Reference 1 along with a description of the calculation procedure. 
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Results : 

Sooting was a serious problem with all of the hydrocarbon fuels tested. 
Methanol did not present this problem and only at very high pressures (around 
1100 psia) did ethanol begin sooting. Propanol- 1, on the other hand, exhibited 
sooting problems at much lower pressures. The upper limit due to sooting for 
n-decane was 32 atmospheres. 

The tests using ethanol as a fuel were terminated around 72 atmospheres 
since it was impossible to determine the nature of the sphere surface at pres- 
sures in excess of this. Methanol presented this problem also but at slightly 
higher pressures. In both cases vapor jets were observed rather than liquid 
drops. Critical combustion could account for this effect although the pressure 
at which this condition occurred could not be precisely determined. 

The experimental burning rates for the alcohols are shown in Figure 3. 

The burning rate data for the hydrocarbon fuels (Figure 4) were terminated 
at high pressure due to the formation of soot. When this occurred, carbon 
spots would form and grow on the surface of the sphere. As can be seen 
from Figures 3 and 4, the theory gives a reasonably good indication of the 
rate of increase of the burning rate with increasing pressure. 

The effect of varying sphere size is examined for methanol, n-heptane 
and n-decane in Figure 5. For this plot, the dimensionless burning rate, 
normalized by the convection correction is employed for the ordinate so 
that the data for various sphere sizes should fall on a single curve. It 
is seen that the normalized burning rate (which corresponds to no-convection 
burning rate of the theory) is almost a constant up to the critical burning 
condition for the present porous sphere experiments. This indicates that 
the burning rate increase with increasing pressure in Figure 2 and 3 is 
largely due to convection effects. The present experimental results re- 
present a reasonably good test of the burning rate correction for natural 
convection, since the Grashoff number based on the Spalding definition (8) 
varies in range 10 4 - 10 8 . 

Liquid surface temperatures for the six fuels are illustrated in 
Figures 6 and 7. The boiling point curves and the surface temperature pre- 
dictions of both the low-pressure and quartenary high-pressure theories are 
shown on the figures along with the data. The data for ethanol and propanol- 1 
(Figure 6) agrees reasonably well with the high pressure theory at the high 
pressures. However, the low pressure theory gives the best estimation of 
data over the entire test range for methanol . The poorer high-pressure 
theoretical results for methanol could be attributed to the large quantities 
of water vapor in the combustion products of this fuel. Water is difficult 
to model precisely in the high-pressure phase equilibrium analysis, and 
materials with high water vapor concentrations have generally shown poorer 
agreement with the high-pressure theory in the past (15), 

The data for the paraffins (Figure 7) could not be extended to suf- 
ficiently high pressures to provide an adequate test of high pressure 
theory due to the formation of soot. However, the low pressure theory 
appears to be adequate for these fuels over the available experimental 
range . 

The simplified binary model gave essentially the same results with 
regard to burning rates and liquid surface temperature as the quaternary 
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phase equilibrium model. However, in contrast to high pressure droplet 
combustion where the binary model also gave a good estimation of critical 
burning conditions there were significant differences between the critical 
porous sphere combustion pressures predicted by the two high-pressure theories. 
The critical combustion conditions for all three theories are compared with 
pure fuel critical properties in Table I. In agreement with the experimental 
findings, both the low and high pressure quaternary theories predict critical 
burning pressures on the order of 100 atm for methanol and ethanol. The 
theoretical indication that ethanol should experience critical burning at 
pressures somewhat below methanol is also in qualitative agreement with the 
fact that experimental difficulties in determining burning rates were 
encountered at somewhat lower pressure for ethanol, c.f. Figure 3. 

The theoretical results were obtained with the variable property-variable 
specific heat gas phase analysis, using the properties listed in Reference 1. 
The use of the variable property- constant specific heat and constant property 
gas phase analysis gave essentially the same results, when the respective 
constant properties in each of these cases was evalutaed at average conditions 
in each region. 


Summary and Conclusions: 

The low and high pressure theories gave essentially the same pre- 
diction of burning rates at high pressures. The greatest differences 
between the theories were encountered in the prediction of liquid sur- 
face temperatures and critical burning conditions. The quaternary high 
pressure theory gave the best agreement with experimental results for 
ethanol and propanol- 1, whereas the low pressure theory gave the best pre- 
dictions for methanol. 

The experiments indicated that the methanol and ethanol were approach- 
ing critical combustion conditions at pressures on the order of 80-100 atm. 
Both the low pressure and high pressure quaternary theories predicted 
critical burning in reasonable agreement with these results. Critical 
burning conditions could not be approached for the remaining fuels due to 
the formation of soot deposits on the sphere over the pressure range of 
30-60 atm. 

The use of variable property- variable specific heat, variable pro- 
perty-constant specific heat and constant property gas phase analyses 
gave essentially the same results as long as the respective constant 
properties were evaluated at average conditions inside and outside the 
flame. Parametric property variations caused variations in the computed 
results similar to those encountered in earlier high pressure combustion 
studies (3, 4, 15). For porous spheres , the binary high-pressure theory 
gave a poorer approximation of the quaternary high-pressure theory, then 
was the case for high pressure droplet combustion (3). 

Work has been completed on this phase of the investigation for droplet 
combustion in a cold gas environment. 
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Nomenclature, Section II: 



c 


C . 
Pi 


D 

g 


K, K 


n. 


n 


o 


P 

r 

R 

e 

T 

X 

v 

X 


Specific heat parameters. Equation (1) 

Concentration 

Specific heat of species i 

Droplet diameter 
Diffusion coefficient 
Acceleration of gravity 
Binary interaction parameters 

Burning rate constant with and without convection 

Total molar flow rate with and without convection 

Prandtl number 

Reynolds number 

Temperature 
Thermal conductivity 
Kinematic viscosity 
Parameter, Equation (3) 


Subscripts 

o Initial conditions 

s Liquid surface 


III. High Pressure Droplet Burning in Combustion Gases 


Introduction: 

The objectives of this portion of the investigation include the 
determination of high pressure droplet burning and evaporation rates 
and liquid surface temperatures in a combustion chamber environment. 
The range of fuels to be considered will be the same as the high pres- 
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sure natural convection tests. The oxygen concentration in the combustion 
gases will be varied and will include the limiting case of high temperature 
evaporation without combustion. In addition, the experiments will include 
various flow velocities past the test droplet and various pressure levels 
to provide an evaluation of both forced convection and high pressure effects. 
The experimental results will also be compared with high pressure combustion 
theories . 

In order to simulate a combustion chamber environment for the porous 
fuel spheres, a burner is used to produce the hot combustion gases. The 
burner provides an added advantage in that the composition of the gases 
can be accurately controlled. The burner design is of the swirl-stabilized 
diffusion type. Other major advantages of a burner of this type include a 
relatively simple design, stable operation over a wide range of flow con- 
ditions, and the elimination of flashback. In earlier tests with premixed 
flat flame burners, flashback posed serious problems. 

The burner used in the current tests employs a flame supplied with 
a mixture of carbon monoxide and air. Carbon monoxide was chosen for a 
fuel since it is quite stable, has a reduced tendency to smoke and minimizes 
the problem of water condensation within the apparatus. 


Experimental Apparatus : 

A sketch of the high pressure combustion rig with the burner is shown 
in Figure 8. The basic reactor vessel is similar to the one that was used for 
the cold gas droplet burning tests. There are, however, several modifications 
that were not present on the former apparatus. These include an internal 
stainless steel cooling coil to keep the apparatus' walls cool and to reduce 
radiation from the walls to the fuel droplet. Observation of the combustion 
process is made through quartz windows located at the test section. The 
fuel delivery and measuring system is the same as that described for the 
cold gas tests. There are four major burner components. These include the 
stainless steel burner base containing the fuel and oxidizer injectors, a 
spark -ignition system, a ceramic flow straightener, and a steel outer jacket. 

The internal diameter of the burner test section is 5.08 cm. The 
burner wall consists of alundum tubing and alumina firebrick. Fuel is 
admitted through an injector located in the center of the burner base 
and the oxidizer is admitted through two tangential injectors located 
along the inner wall of the burner. A flow straightener is located in the 
burner passage to insure that the velocity and temperature profiles of the 
hot gases will be uniform when the test section is reached. All gas flows 
are metered by critical flow orifices. To determine gas temperature, a 
system of platinum-platinum rhodium fine wire thermocouples is located in 
the test section. To reduce oxidation and catalytic effects the thermocouples 
are coated with silicone using a procedure outlined by Fristrom and Westen- 
berg (11). As a check of the gas temperature measurements, the composition 
and temperature of the gas at the droplet location are also computed from 
thermochemical calculations. For these calculations, allowance is made for 
dissociation and heat loss through the burner walls. 
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Theoretical Analysis: 

The theoretical results and analysis of the high pressure theory will 
be employed in this phase of the investigation. As was pointed out in the 
discussion of results in the portion of the report dealing with high pres- 
sure droplet combustion, the high pressure model gave adequate agreement 
with the experimental results for the range of fuels tested. Therefore, 
this model will be used but with one major exception. The convection cor- 
rection that accounts for the influence of the flow of burning gas past the 
test droplet is different than that used for the cold gas tests (which were 
essentially free convection tests) . 

The specific correlation for the present tests follows a fojm which 
was previously employed in Reference 5. 

K/K q = 1 + 0.278 R e l/2 Pr l/s (1 + 1 . 237/RePr 4/3 ) ~ l/l 

where K and K q are the burning rate constants with and without convection present . 

This expression has the advantage of asymptotically approaching the 
mathematically rigorous results of Fendell et.al (12) for small Peclet 
numbers as well as agreeing with the data of Frossling (13), Yuge (7) and 
Allender (14) for 10<Re<800. 


Results : 

The experimental apparatus for high pressure droplet burning rate 
measurements in a combustion gas environment has been assembled. In 
preliminary tests, the swirl stabilized burner has been operated for 
extended periods of time at elevated pressures. During these tests, no 
difficulties were encountered with regard to overheating the structure 
of the apparatus. , The burner flame was found to be quite stable for oper- 
ation at various gas pressures, mixture ratios and velocities. 

Testing with porous spheres is currently in progress. Complete 
combustion testing as well as analysis of the results will be completed 
during the coming report period. The pressure range of these tests will 
extend to pressures on the order of 35 atm, which is sufficient to reach 
conditions where critical point phenomena become important. 


Summary and Conclusions: 

An apparatus has been constructed to allow measurements of steady 
evaporation and combustion of liquid fuels in an environment which 
simulates combustion chamber conditions. During the coming report period, 
testing will be completed for a range of fuels, oxygen concentrations, 
gas velocities and pressures; these results will be compared with analytical 
results which allow for the effect of high pressure phenomena, ambient gas 
solubility, etc. 


\ 
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IV. Liquid Strand Combustion 


Introduction: 

The purpose of this portion of the investigation is to study the temperature 
distribution and dissolved gas concentration in a burning liquid monopropellant 
strand at high pressures. Earlier work concentrated on the nitrate ester fuels 
(15). The study described here is concerned with hydrazine. 

During this report period, burning rate measurements of liquid strands 
of hydrazine were extended from previous by reported results (16) to include 
the effect of tube diameter and liquid phase concentration. Also, liquid 
surface temperature measurements were obtained over a pressure range of 7 
to 300 psia. These temperature measurements were compared to the infinite 
activation energy theoretical model of Faeth (15) and also to predictions 
based on the familiar unity Lewis number assumption. 


Experimental Apparatus : 

Since the strand burner apparatus has been described in detail in 
Reference 16, only a brief description of its operation will be given here. 
With the strand burner technique, liquid fuel is placed in a glass tube con- 
tained within a windowed pressure vessel. The vessel is pressurized with 
nitrogen to the desired test pressure. Following ignition with a heater 
coil, the fuel bums down the tube. As the combustion zone propagates past 
the window in the chamber, the rate of regression of the liquid column 
as well as the position of the thermocouple with respect to the liquid 
surface is determined from motion picture shadowgraphs. The temperature 
registered by the thermocouple is recorded on an oscillograph with a flat 
frequency response to 2000 Hz. 

The thermocouples were constructed of 0.0003 inch O.D. platinum - 
platinum 10% rhodium butt welded wire as described in Reference 11. Since 
the pressure levels and burning rates for the hydrazine tests were consider- 
ably lower than the conditions encountered in earlier testing (15, 16), the 
cold reference junction of the thermocouple was relocated outside the ap- 
paratus in order to obtain a steady reference temperature. Otherwise, the 
apparatus used was as described in Reference 15, 


Theoretical Analysis: 

Since the process is steady, the gas and liquid phases are considered 
separately. However, phase equilibrium requires the two solutions to match 
at the liquid surface. 

Two gas phase models were used. The first model has been described in 
detail by Faeth (15), This model assumes an infinitely thin decomposition 
flame (infinite activation energy) and the gas phase specific heat and 
thermal conductivity are assumed to vary linearly with temperature. In the 
second model, the unity Lewis number assumption is used and the gas phase 
specific heat is assumed to be constant but thermal conductivity is as- 
sumed to vary linearly with temperature. For both models the ratio of gas 
phase thermal conductivity to the concentration of fuel times the diffusion 
coefficient is assumed to be constant. The flame temperature, assuming an 
adiabatic deflagration wave, and product gas concentrations were determined 
from thermochemical calculations allowing for all relevant dissociation 
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reactions. Data for these calculations were obtained from the JANNAF Tables (17). 
All other properties were computed as in Reference 1. 

Two models of phase equilibrium at the liquid surface have previously 
been used (15) . The low pressure model neglects high pressure corrections 
and ambient gas solubility in the liquid phase. The mole fraction of fuel 
at the liquid surface is taken to be the ratio of vapor pressure of fuel 
at the liquid surface temperature to the total pressure. The more complete 
high pressure model accounts for high pressure effects and the solubility of 
the combustion products in the liquid phase. Since the pressure levels for 
the hydrazine tests were much less than the critical pressure, only the low 
pressure phase equilibrium model has been applied to the hydrazine data. 

With the low pressure phase equilibrium model, one relation between the 
liquid surface temperature and fuel mole fraction at the surface was 
obtained from the vapor pressure curve for a given total pressure. A 
second relation was obtained from the gas phase model (either infinite 
activation energy or unity Lewis number) . The point where the two relations 
were simultaneously satisfied gave the unique solution for the liquid surface 
temperature and surface fuel mole fraction for a given total pressure. 


Results : 

The burning rate results for hydrazine are summarized in Figure 9. For 
the highest purity hydrazine tested, the burning rate varied linearly with 
pressure for all tube sizes except for the lowest pressures tested where 
a slight fall off from the linear relation occurred. As pressure increased 
the tube size had less effect on the burning rate but the liquid concentration 
had a significant effect. 

The purity of the hydrazine sample was obtained by commercial gas 
chromatographic analysis. This analysis indicated that the sample was 
made up of 98.6% hydrazine, 1.3% water and 0.1% trace impurities, by 
volume. The purity of the other samples of hydrazine shown on Figure 9 
was obtained by water dilution of the 98.6% sample. 

The upper limit in pressure shown on Figure 9 was the highest pressure 
for which hydrazine would steadily propagate. Upper extinction pressures 
for the stand burning of hydrazine have previously been noted by other 
investigators (18, 19, 20). As indicated on Figure 9, the upper extinction 
pressure varied with tube size and also liquid concentration. 

The liquid surface temperature results are shown on Figure 10. As 
expected, the low pressure phase equilibrium model was quite adequate for 
hydrazine since the pressure levels tested were well below the critical 
pressure. Also evident from Figure 10 is that the two different gas 
phase models predict nearly identical results and both predictions agree 
quite well with the experimental values. 


Summary and Conclusions: 

The burning rates and surface temperatures of liquid strands of hy- 
drazine have been measured for various pressures. Burning rates were 
obtained over the pressure range of 7 to 600 psia and were found to in- 
crease linearly with pressure except for subatmospheric pressures. The 
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concentration water in liquid hydrazine has a significant effect on the burning 
rate particularly at high pressures. Smaller tube diameters were found to 
cause an increase in the apparent burning rate, particularly at low pressures. 

Surface temperatures of burning liquid strands of hydrazine were measured 
in the pressure range of 7 to 300 psia. The low pressure phase equilibrium 
model with and without using the unity Lewis number assumption predicts the 
experimental data quite well over the entire pressure range. 


V. Oscillatory Combustion 


Introduction : 

This portion of the investigation concerns the influence of pressure 
transients on the combustion of a liquid monopropellant fuel. The com- 
bustion system consists of a liquid surface from which fuel evaporates; 
as the gaseous fuel flows away from the liquid surface, it undergoes 
exothermic chemical decomposition. A portion of the heat evolved in the 
decomposition process flows back to the liquid surface providing the 
energy for further evaporation of the liquid. The remainder of the 
energy released is convected downstream by the gas flow. The main thrust 
of the present work is to gain an understanding of the coupling between 
an oscillatory pressure wave and the fuel burning rate. 

Numerous theories have been developed to predict the response of a 
burning liquid to pressure oscillations. However, direct experimental 
measurements of this phenomenon are lacking. The purpose of this study 
is to obtain experimental data on the response of a burning monopropellant 
to an oscillating pressure environment and to compare these results with 
a theoretical model of the system. 

An important criterion must be met if an adequate test of any theoret- 
ical response model is to be made. For a burning liquid, the characteristic 
frequency of thermal wave propagation in the liquid phase can be defined 
as m 2 /a where m is the mean burning rate of the fuel and a is the liquid 
phase thermal diffusivity. For pressure oscillation frequencies much less 
than this characteristic frequency, the system behaves quasi-steadily . 

This frequency regime provides a relatively weak test of theoretical response 
models. Therefore, pressure oscillation frequencies should be at least the 
same order of magnitude as the liquid phase critical frequency for at least 
a portion of the experimental data. 

The steady strand study indicated that the burning rate, m, of 
hydrazine is roughly proportional to pressure. Since a for hydrazine is 
nearly independent of pressure, the characteristic frequency decreases 
with decreasing pressure. Quasi-steady analysis indicated that low mean 
pressures produce larger surface oscillations for a given pressure 
amplitude (16) . Thus low pressures are favored as the condition where the 
response of the liquid can be observed throughout the entire frequency 
range of interest. Since hydrazine will burn as a liquid strand at low 
pressures, this fuel is particularly suited for measurements of burning 
rate response and it was chosen as the fuel to be studied. 
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During this report period an oscillatory pressure version of the 
steady strand apparatus was assembled and preliminary testing was begun. 
Also, the theoretical model was formulated and the resulting equations 
derived. These equations are now being coded for computer solution. 


Experimental Apparatus : 

The oscillatory combustion apparatus is shown in Figure 11. The major 
components are the glass tube filled with liquid monopropellant fuel, 
the camera- lens system used to measure the fluctuations of the burning 
liquid surface to the imposed pressure oscillations and the rotary valve 
arrangement used to provide the oscillatory pressure variations. The fuel 
ignition system is similar to the one used in the steady strand study. 

The rotary valve arrangement is similar to the system described in 
Reference 21 for establishing an oscillating propane gas flame. With this 
technique an oscillating pressure is set up in the test chamber by an 
oscillating air stream. 

The oscillating air stream is produced through the use of a ball 
valve and needle valve mounted in parallel. The amplitude of the oscillat- 
ing pressure can be varied independently of the frequency by varying the 
relative amounts of air passing through the two valves. The frequency is 
varied by varying the speed of the DC motor used to rotate the ball valve. 
The mean pressure in the chamber can be varied by the setting of the needle 
valve in the exhaust line. For testing at subatmospheric pressures, the 
exhaust line is connected to a large chamber which is in turn connected to 
a vacuum pump. 

Quasi-steady analysis indicates that typical liquid surface oscilla- 
tion amplitudes are on the order of 0.01 mm (16). In order to resolve such 
small distances, the camera-lens system shown on Figure 11 is qsed in con- 
junction with a Vanguard motion analyzer. The camera is a 35 mm strip 
camera and the lens has a focal length of 150 mm. Primary gains as high 
as 5:1 have been obtained with this system. This gain coupled with the 
gain of about 20:1 of the Vanguard provides an overall gain of 100:1 
which is adequate for accurate resolution of the film records. 

Any phase shift between the imposed pressure oscillations and the 
resulting burning rate oscillations are determined from the films and the 
pressure trace registered on an oscillograph. The film and pressure re- 
cords are synchronized by a switch closure which deflects a galvanometer 
on the pressure record and starts a light streak on the film. Timing 
marks are placed on both the film and pressure record. 


Theoretical Analysis: 

Existing Theories 

A number of theories have been developed to determine the res- 
ponse of burning liquids to imposed pressure and velocity fluctuations 
(22-24, to name only a few). However, none of the current liquid sur- 
face response models provide for a monopropellant gas phase combustion 
process. Earlier work on the combustion of hydrazine indicates that 
a monopropellant gas phase model is more realistic for the hydrazine 
fuels under rocket engine conditions (25) . 
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Since the process is one dimensional, some insight into the 
modeling of the experiment can be gained by studing solid propel- 
lant combustion instability theories (26-31 are representative). 

In particular, neglecting condensed phase reactions, solid propel- 
lant models are applicable to the present problem with modified 
boundary conditions at the two phase interface and slight modifica- 
tions of the gas phase analysis. Although the present theoretical 
model is of the one-dimensional experiment, the results should lead 
to a better understanding of liquid response and thus to improved 
droplet response models. 


General Model 


A sketch of the theoretical model is shown on Figure 12. The 
coordinate system chosen is inertial with the origin at the mean 
position of the burning surface (28-29) . With the imposed sinusoidal 
pressure oscillations, x , the instantaneous position of the liquid 
surface, varies as shown in Figure 12 to first order in pressure ampli- 
tude. The steady state temperature, fuel mass fraction and gas phase 
reaction rate profiles are also shown on Figure 12. 

For generality, the gas phase transient effects have been included 
in the analysis similar to the approach used by T'ien for solid pro- 
pellants (31) . Also included are the effects of variable properties 
through the introduction of a modified Howarth transformation. 

The major assumptions used in the analysis are as follows: 

1. The flow is one-dimensional with a Mach number much less 
than unity and all body forces are negligible. 

2. The flame is premixed and laminar. A one-step, irreversible 
chemical reaction takes place in the reaction zone and any 
time lags associated with the chemical reaction itself are 
negligible . 

3. Radiation heat transfer is negligible. An estimation of the 
radiation contribution indicates that it accounts for less 
than 1% of the total energy required for the vaporization 

of hydrazine under steady state conditions. This is not 
surprising since hydrazine has a low adiabatic flame tempera- 
ture . 

4. With regard to thermodynamic properties, it is assumed that 
all gas phase diffusion coefficients are equal, all molecular 
weights are equal and constant, all gas phase specific heats 
are equal and constant, the gas phase mixture thermal conduc- 
tivity varies linearly with temperature but is independent of 
composition, and the liquid is composed of a single chemical 
species with constant properties. 

5. Since the pressure levels are low, the combustion products 
are taken to be insoluble in the liquid phase. Chemical 
reaction is also neglected in the liquid phase. 

6. Based on the steady strand study, the Lewis number is assumed 
to be unity. 
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7. The gas phase is taken to be an ideal gas and at the liquid 
surface the surface fuel mass fraction is related to the 
surface temperature through the Clausius-Clapeyron equation. 

8. Since the Mach number is much less than unity, the inertial 
and viscous terms in the momentum equation are neglected. 

9. The complete unsteady gas phase equations are used but the 
wavelength of any periodic pressure disturbance is assumed 
to be long compared to the distance from the outer edge of 
the reaction zone to the liquid surface. Estimation 
indicates that this assumption is valid for frequencies up 
to about 50,000 Hz. Assumptions 8 and 9 indicate that pres- 
sure is only a function of time. 

Thus, the theoretical model includes the effects of time depend- 
ent heat and mass transfer in both the liquid and gas phases. Vari- 
able property effects are also included in the gas phase. Global 
kinetic parameters are used to characterize the hydrazine decomposi- 
tion flame. 


Governing Equations 

With the above assumptions, the equations of overall continuity, 
species conservation and energy in the gas phase and overall contin- 
uity and energy in the liquid phase were obtained from the general 
equations given by Williams (32). 

For the gas phase the equations are* : 


Gas Phase, x g *(t*) < x* < °° 


Conservation of Mass 


8p* 3p*v* 

3t* + 


3x* 


CD 


Conservation of Species 


3Yi 

P* at* + 


p*v* 


3Y . 

i_ 

ax* 


a 

ax* 


p*D 


3Y. 

l 

ax* 


'.w* c 
i F 


( 2 ) 


Conservation of Energy 


p*C. 


8T* 


+ p*v*C 


3T* 


P 9t * P 9x* 


a 

ax* 


3T* 

ax* 


ap* 

at* 


q*w* 


(3) 


Since the gas is ideal 


p* = p*RT* 


( 4 ) 


*See the nomenclature at the end of this section for the meaning 
of symbols. 
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In Equations (2) and (3) the following definitions have been used: 


w 


* 

F 


q* 


/. 

l 


B*T< 


p* ' 

n 

[V 

n 

r \ 

-E* 

RT*j 


Im fJ 

exp 

[rt*J 


n /. M. 

-E h°-^ 

i-1 1 /. M„ 

l F 

i 


(4 


/ 


A) M'. 
1 1 

M„ 


(5) 

( 6 ) 


(7) 


•k 

Here w p is the^rate of reaction of the fuel with a pre-exponential 
factor of B*T* and an activation energy of E* , q* is the heat of 
combustion and /. is related to the stoichiometric coefficients. 

The equations of. conservation of mass and energy in the liquid 
phase are: 


Liquid Phase, -« < x* < x g * (t* ) 


Conservation of Mass 


v^* * constant 


( 8 ) 


Conservation of Energy 


p*C 


El 

£ 8t v 


* * 3 j* 

P £ V £ C £ Jx* 


- X, 


8 2 T* 

* r 

Sx * 


= 0 


(9) 


Since only oscillatory solutions are required, the initial condi- 
tions are irrelevant. The boundary conditions applicable to the pre- 
sent problem can be summarized as follows. At the cold end of the 
liquid propellant, the temperature must be constant with respect to 
both x* and t* . 


X* -V -00 

T* "*■ T* (a constant) (10) 

At the liquid surface x *(t*) six conditions must be satisfied. The 
first two are obvious and are mentioned only for completeness. That is, 
the temperature and pressure must be continuous across the surface. 


T* 

s 


= T* 


P \ 

s 


( 11 ) 


( 12 ) 
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Conservation of riiass applied across the moving surface reduces to: 


* * 


P„ Cv„ - = p*(v* - x s ) 


( 13 ) 


and conservation of energy at the surface yields: 


3T*i 


* 8T*' 

9x*j 

s 

X 91c" 


- «X - v its - c iK * L * ] C143 


Since the products of combustion are assumed to be insoluble in the 
liquid phase, the gradient of the fuel mass fraction is related to 
the fuel mass fraction at the liquid surface as follows: 


3Y F 
P * D 9lF 


★ * 
‘ > t Cv * 


V ' V F,s - » 


(15) 


The Clausius-Clapeyron equation relates the fuel mass fraction at the 
surface to the surface temperature. 



^ exp 


-L* 


RT 


(16) 


where L* is the heat of vaporization of the fuel and a* is a constant 
with the dimensions of pressure. The remaining boundary conditions re- 
late the variation of fuel mass fraction and temperature far from the 
liquid surface. Far from the surface, the gas phase chemical reaction 
must go to completion since eventually all the fuel will react. Thus, 


X* OO 

V p - 0 (17) 

for all time. In addition, the temperature must become independent 
of x* since the only energy source, the chemical reaction, has gone 
to completion and heat conduction will smooth out any temperature 
variations with respect to distance. . However, the temperature will 
still vary with time since the pressure is a function of time. There- 
fore : 


x* <» 

T* (x* , t*) + T* (t*) (18) 

where T*(t*) is a known function. 

oo v J 

Equations 1-18 are placed in dimensionless form by introducing 
the following variables. 
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T 


T* 

-p* 

OOQ 


V = 



P 


P* 

P* 

o 


A 


A* 

OOQ 


t 


t*p* C v> 


OOQ 


00 o 


A* 

«o 


(19) 


x 


x*p* C v* 

ooo p 000 

A* - 


OOQ 


The subscript °°o signifies a zero order (steady state) quantity evaluated 
at infinity. 


Treatment of variable gas phase properties is accomplished by the 
Howarth Transformation appropriate to the present coordinate system by 
defining 


n 


C v* 


OOQ 


A* 

°°o 


p*dx* 


x*(t*) 

s 


( 20 ) 


and by defining a modified mass flux as 


r 


pv + 


3_ 

at 


(X 

pdx 


x s Ct) 


( 21 ) 


The Shvab-Zeldovich variable 9, is also introduced as a dependent 
variable , 


6 - qYp + T (22) 

Substituting equations 19-22 into equations 1-4 then yields: 


Gas Phase, 0 < n < 00 


3r 

an 

0 





(23) 

ae 

ae 

p 329 

Y-l 

1_ 

3P 

(24) 

— + 
at 

r an 

an 2 

y J 

p 

at = 0 

3T 

8T 

p ill 

Y-l'' 

1_ 

ap 

(25) 

- ■■ + 

at 

r an 

an 2 

V 1 

p 

aT = qw 

p = 

pT 





(26) 


In deriving Equations (24) and 
that /p = -1 for the present case. 


(25) use has been made of the fact 
The dimensionless reaction rate is 
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defined as 


A_ x 1 + 6 

q 



(6 - T) n exp 



(27) 


where E is the dimensionless activation energy. 


E 


RT* 


CO o 


(28) 


q is the dimensionless heat of combustion 


q (29) 

p °°o 

and A is related to the pre-exponential factor and zero order mass 
flux 


A 


B*T* A* 

»0 °°0 

p 2 v 2 C 
°°o °°o p 


M„RT* 

F °°o 


n 


(30) 


Also, noting that 


Af _ T* 

A* T* 

ooo °°0 


(31) 


and through the use of the ideal gas equation of state, it follows that 

pA = P (32) 

The dimensionless liquid phase continuity and energy equations are 
Liquid Phase, -<*> < n < 0 

r = r(t) (33) 


3T 3T 

3t + r 9n 



= 0 


The dimensionless parameter 6^ 

6, 


P £ X £ C p 


is a constant defined as 


(34) 


(35) 


The nondimensional boundary conditions are as follows. At the cold 
end of the fuel, 

n -> -« 

T(n,t) - T_ m (36) 
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At the liquid surface, the conservation of mass yields 

n = 0 


r(t) = p (o) [v(o)-x s ] = P £ (v £ -x s ) 


( 37 ) 


The conservation of energy across the surface becomes 

91 ] 


P lh 


3T 

3n 


n = o 


( 3r U 


n = o 


+ - r[T(o) (1-g) + L] 


where 3 is the ratio of liquid to gas specific heats 

6 - ^ 

6 " c 


(38) 


(39) 


and L is the nondimens ional heat of vaporization 
I * 

I = L 
C T* 
p °°o 

The insolubility condition requires that 


30 


+ = r [0 (o) - T (o) - q] + P 


n = o 


3T 


3n 


n = o 


(40) 


(41) 


The Clausius-Clapeyron equation after nondimensionalization is 


0(0) = exp + T (°) 

where 

a* 

a - P * 
o 

The dimensionless conditions far from the surface reduce to 


(42) 


(43) 


n -*■ °° 

e(n, t) T(n,t) + T(t) (44) 

where T(t) is a specified function. 

To solve the system of partial differential equations, Equations 
(23-25) and (34), with the given boundary conditions, it is assumed 
that all amplitudes of oscillation are small, allowing the use of a 
perturbation analysis. For oscillatory solutions, the dependent 
variables can be expressed as 
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P(t) = 1 + ee la)t 

T(n, t) = T o (n) + eT 1 (n) e lwt 

e(n, t) = Q 0 (n) + ee 1 (n)e iwt (45) 

r(t) = r + er-e^ 1 " 

v J o 1 

x (t) = 0 + ex 1 e : '' a3t 

s si 

where e represents the normalized amplitude of the imposed pressure oscil- 
lations. Substituting Equation (45) into Equations (23-25) and (34) and 
the boundary conditions, the resulting equations are separated in terms 
of powers of ee lait . 


Zero Order (Steady State) Problem 

For the zero order problem, the gas phase can be solved independ- 
ently of the liquid phase. The problem reduces to the following: 

Gas Phase, 0 < n < 00 


^2^ clT 

° + q l-" AT 1+6 " n (1-T ) n exp 

dn 2 dn ° ° 


-E 


= 0 


(46) 


with the boundary conditions 
n -*■ °° 


T -► 1 
o 


and at n = 0 


(47) 


fdT 1 
o 


dn 


n = 0 


* = q + T o (D) - 1 


-L 


T (0) = 1 - aq exp [ ? fQ J 


(48) 

(49) 


In Equations (46-49) use has been made of the conservation of mass require- 
ment that 


r = constant = p v = 1 

O 00 O OOQ 


(50) 


and also that 


e 0 M - l 


( 51 ) 


from Equations (24) and (45) . 
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In Equation (46) , A is unknown since it is related to the burning rate 
p* v* o which is unknown. The three boundary conditions permit a unique 
va?ue of A to be determined. 


First Order Problem 

The liquid phase solution is required for the complete specification 
of the first order problem. To zero order in e, the equation and bound- 
ary conditions are: 


Liquid Phase, -» < n < 0 


6 £2To _ dTo 
^dn 2 dn 


0 


n -*■ -°° 


T 

o 


-* T 

_ CO 


n = 0 


T 

o 


T o (°) 


(52) 


(53) 


The solution of Equation (52) subject to the boundary conditions of 
Equation (53) is 


r o = T .» * t T 0 C°) - T _J e*P Cn/«P 


(54) 


To first order in e the equation and boundary conditions for the liquid 
phase are: 


- dfl 

dn 2 


dn 


iuT, = r 


dT 

c 

1 dn 


(55) 


n -*■ -°° t, o 

n = 0 

The solution of the first order system is 


T x = T x (0) 


T 1 ■ {T 1 (0) * 6^ lT o (0) - T -J> 


r 1 + /1+4. 6, 

ICO 1 


26. 


(56) 


[T 0 COJ - T_J exp (n/sp 


(57) 
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The first order gas phase equations reduce to the following. 
Gas Phase, 0 < n < 00 


r^ = constant 


( 58 ) 
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In Equation (60) w q is defined by the following expression: 


A „ 1 + 6 - n „ .n 

w = — T (1-T) exp 

o n o o * 

q 


-E 
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(61) 


Using Equation (57) the conservation of energy boundary condition 
to first order in e becomes 
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The Clausius - Clapeyron equation reduces to 
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The outer boundary conditions from Equation (47) are 

n 00 

0 j Cn) T^(n) -> a constant (65) 

The constant in Equation (65) can be determined by solving Equation 
(59) . It is easily verified that as n -»■ 00 , 


6 1 (n) Tj (n) -> (66) 

Y 

which is the form for isentropic flow. 

Now consider how the theoretical results can be compared to the experiment 
results. For the zero order (steady state) problem, the steady strand 
measurements provide steady state burning rates and liquid surface temperatures. 
By matching theory with experiment, realistic global kinetic parameters 
(activation energy and pre-exponential factor) can be determined for given 
thermodynamic properties. 

For the unsteady case, e, the expansion parameter, is related to the 
amplitude of the imposed pressure oscillations. Since r^ is related to the 
unsteady burning rate in response to these imposed pressure oscillations, 
it remains to be shown how r^ is related to the unsteady burning rate 
measured in the experiment. 

From the conservation of mass at the liqud surface. Equation (37), and 
Equation (45) it follows that 

r l - -<’* 1 ™‘sl (67) 


Rearranging Equation (67) 

- r l 
x 


si 




( 68 ) 


Since is a given thermodynamic property, x ^ is directly related 
to r^. Moreover, x ^ is a complex number with botfi amplitude and phase 
parameters. Both of these parameters can be directly compared with the 
measured amplitude of the liquid surface oscillations and the phase lag 
of these oscillation with respect to the pressure oscillations. 


At the present time the equations are being coded for computer 
solution. The zero order theoretical solution will be matched to the 
steady strand work to determine realistic global kinetic parameters. 

Using these same parameters, the first order equations will be solved to 
determine the predicted burning rate f luxuations . These predicted values 
will then be directly compared to the unsteady experimental burning rate 
results . 
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Summary and Conclusions: 

The experimental apparatus has been assembled and preliminary testing 
has begun. During the coming report period, the burning rate response of 
hydrazine will be measured for various imposed pressure wave amplitudes 
and frequencies of oscillation. The frequency range tested will include, 
at least, frequencies of the same order of magnitude as the liquid phase 
critical frequency. In nondimensional form this corresponds to u on the 
order of 10 ^ to 10 2 . 

The oscillatory theory has been developed through first order in 
amplitude of the pressure oscillations. The resulting equations are being 
coded for computer solution. The theoretical predictions will be compared 
to the experimental results to provide a direct test of the model . 


Symbol 

A 


a 

B* 



i 


L 

M 


m 

N 


n 


Nomenclature, Section V 
Description 
Defined in Equation (30) 

Constant in Claus ius-Clapeyron Equation, see Equation (16) 

Pre-exponential factor 

Specific heat of liquid 

Specific heat of gas 

Diffusion coefficient 

Activation energy 

Standard heat of formation of species i 

/T 

Heat of vaporization 
Molecular weight 
Mean burning rate 
Total number of species 
Reaction order 



Symbol 

P 

q 

R 

r 

T 

T 

— CO 

t 
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w 

X 

x s 
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a 
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Y 
6 

6 1 

e 

n 
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Description 

Pressure 

Heat of combustion 
Gas constnat 
Defined in Eq. (21) 

Temperature 

Temperature of cold end of liquid fuel 

Time 
Velocity 
Reaction rate 
Distance 

Instantaneous position of liquid surface 
Mass fraction 
Liquid diffusivity 

Ratio of specific heats, see Eq. (39) 

Specific heat ratio 

Temperature dependence of pre-exponential factor 
Defined by Eq, (35) 

Amplitude of oscillatory pressure 
Dimensionless distance, see Eq. (20) 

Defined by Eq. (22) 

Thermal conductivity 
Defined by Eq. (7) 

Stoichiometric coefficient of species i as reactant 
Stoichiometric coefficient of species i as product 
Density 

Nondimens ional frequency 
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Subscripts 

F 

i 

a 

s 

s" 

+ 

s 

0 

1 

oo 


Supercripts 

* 


Description 

Fuel 

Species i 
Liquid 

Liquid surface 

Inner side of liquid surface 
Outer side of liquid surface 
Steady State quantities 
First order quantities 
Far from liquid surface 


Dimensional quantities 
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Predicted Critical Burning Conditions for Porous Sphere Combustion in Air 



Fuel inlet and ambient air temperature of 300°K, pressure in atm, temperatures in 
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FIG. 1 HIGH PRESSURE COMBUSTION APPARATUS. 
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FIG. 2 MODEL OF THE BURNING OF A FUEL DROPLET IN AN OXIDIZING 
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